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Abstract 


We  show  how  to  use  polynomial  and  strongly  polynomial  capacity  seeding  algorithms  for 
the  transshipment  problem  to  design  a  polynomial  dual  network  simplex  pivot  rule.  Our 
best  pivoting  strategy  leads  to  an  0(m^  log  n)  bound  on  the  number  of  pivots,  where  n  and 
m  denotes  the  number  of  nodes  and  arcs  in  the  input  network.  If  the  demands  are  integral 
and  at  most  B,  we  also  give  an  C'(m(m  + n  log  n)min(logn5,m  log  n))-time  implementation 
of  a  strategy  that  requires  somewhat  more  pivots. 
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1      Introduction 


The  transshipment  problem  is  one  of  the  central  problems  in  network  optimization.  There 
exist  several  polynomial  and  strongly  polynomial  algorithms  for  solving  this  problem  (see  the 
surveys  [7,1])-  Nevertheless,  the  method  of  choice  in  practice  still  seems  to  be  the  network 
simplex  method. 

In  this  paper  we  shall  consider  the  dual  network  simplex  method.  We  show  that  some 
excess  scaling  algorithms  for  the  uncapacitated  transshipment  problem  can  be  used  to  guide 
the  pivot  selection  of  the  dual  network  simplex  algorithm  for  both  the  capacitated  and  the 
uncapacitated  transshipment  problems.  The  resulting  simplex  algorithm  can  be  viewed  as  a 
special  implementation  of  the  capacity  scaling  algorithm.  This  implementation  maintains  the 
property  that  all  arcs  with  non-zero  flow  are  in  a  tree,  a  property  that  seems  to  help  the  practical 
performance. 

We  give  a  simple  pivoting  strategy  that  leads  to  an  0{n^  \ognB)  bound  on  the  number  of 
pivots  for  the  uncapacitated  transshipment  problem,  assuming  that  the  demands  are  integral, 
and  at  most  B.  We  also  show  how  to  modify  this  strategy  to  achieve  a  strongly  polyno- 
mial 0(0"^  log  n)-pivot  algorithm.  We  describe  a  more  complicated  strategy  that  leads  to  an 
0(nm  log  n)  bound  on  the  number  of  pivots.  The  first  two  pivoting  strategies  can  be  imple- 
mented using  Fibonacci  Heaps  to  run  in  0{n{m  -f  n  log  n)  log  nB)  and  0{Ti^{m  +  n  log  n)  log  n) 
time  respectively.  Dual  network  simplex  algorithms  for  the  uncapacitated  transshipment  prob- 
lem can  also  be  used  to  solve  the  capacitated  version  of  the  problem.  Bounds  for  the  resulting 
algorithms  can  be  obtained  by  substituting  0{m)  for  n. 

Earlier  versions  of  this  paper  have  appeared  as  technical  reports  and  in  conference  proceed- 
ings. The  technical  report  of  Orlin  [9]  described  the  first  polynomial  and  strongly  polynomial 
dual  network  simplex  pivoting  strategies.  These  strategies  were  based  on  capacity  scaling  al- 
gorithms. He  also  gave  a  fcist  implementation  of  the  resulting  0(n^  log  n)-pivot  dual  network 
simplex  algorithm.  The  extended  abstract  of  Plotkin  and  Tardos  [12]  presented  an  improved 
0(nm  log  n)  pivoting  strategy.  The  strongly  polynomial  simplex  algorithm  described  by  Orlin 
in  [9]  assumed  a  model  of  computation,  in  which  we  are  allowed  to  use  some  other  operations 
in  addition  to  the  usual  arithmetic  operations  (addition,  multiplication,  and  comparison).  The 
algorithms  presented  in  this  paper  use  additions  and  comparisons  only. 

The  dual  network  simplex  algorithms  described  in  this  paper  are  based  on  polynomial  and 
strongly  polynomial  excess  scaling  algorithms.  The  scaling  algorithms  work  in  iterations,  where 
each  iteration  executes  an  augmentation  between  a  pair  of  nodes.  The  simplex  implementation 
maintains  a  tree  T  that  contains  all  arcs  with  non-zero  flow  values.  The  tree  is  changed  using 
simplex  pivot  steps,  and  all  augmentations  are  done  through  the  arcs  which  are  in  the  current 
tree.  The  pivoting  strategy  that  results  in  the  lowest  number  of  pivots  is  based  on  Orlin 's  [11] 
strongly  polynomial  transshipment  algorithm.  A  direct  translation  of  this  algorithm  into  the 
simplex  framework  is  infeasible.  For  some  pairs  of  nodes  that  the  algorithm  might  choose  for 
augmentation  there  might  be  no  sequence  of  pivot  steps  that  make  an  augmentation  between 
these  two  nodes  possible  in  the  tree.  The  version  of  the  transshipment  algorithm  used  here 
allows  greater  freedom  in  the  choice  of  the  augmentation  done  at  each  iteration,  and  we  show 


that  a  simple  pivoting  strategy  can  be  viewed  sls  implementing  one  of  the  possible  choices.  In 
effect,  we  let  the  simplex  algorithm  choose  among  the  possible  augmentations. 

A  related  important  open  problem  is  whether  there  exists  a  primal  simplex  pivoting  strategy 
for  the  transshipment  problem  that  leads  to  a  polynomial  bound  on  the  number  of  pivots. 
Some  special  cases  of  the  transshipment  problem  are  known  to  be  solvable  by  polynomial 
versions  of  the  primal  simplex  method.  There  are  polynomial  and  strongly  polynomial  primal 
simplex  algorithms  known  for  the  single  source  shortest  path  problem  and  the  assignment 
problem  (see  [3,10]).  Recently,  Goldfarb  and  Hao  [8]  gave  a  pivoting  strategy  for  the  primal 
network  simplex  method  that  solves  the  maximum  flow  problem  in  0{mn)  pivots  (see  also  [6]). 
Tarjan  [13]  developed  the  first  subexponential  primal  simplex  algorithm  for  minimum-cost  flow 
problem.  The  pivoting  strategy  used  in  Tarjan's  algorithm  is  guided  by  a  polynomial  cost 
scaling  algorithm. 

The  paper  consists  of  four  sections.  Section  2  reviews  the  terminology  and  the  dual  network 
simplex  framework.  Section  3  presents  our  simplest  pivoting  strategy  and  describes  a  modi- 
fication of  this  strategy  that  leads  to  a  strongly  polynomial  bound  on  the  number  of  pivots. 
Section  4  presents  an  improved  strategy  that  leads  to  a  better  strongly  polynomial  bound  on 
the  number  of  pivots. 


2      Preliminaries 


In  this  section  we  define  the  transshipment  problem,  review  some  fundamental  facts  about  it, 
and  review  the  dual  network  simplex  framework.  A  network  is  a  directed  graph  G  =  {V,  E). 
We  shall  use  n  and  m  to  denote  the  number  of  nodes  and  the  number  of  arcs  in  this  graph, 
respectively.  To  simphfy  the  bounds,  we  assume  that  m  >  n.  For  notational  convenience,  we 
assume  that  G  has  no  parallel  or  opposite  arcs.  If  there  is  an  arc  from  a  node  n  to  a  node  w, 
this  arc  is  unique  by  the  assumption,  and  we  denote  it  by  {v,w). 

The  input  to  the  transshipment  problem  consists  of  a  network  G  =  (V,  E),  a  capacity  function 
u  :  E  -^  R"*"  U  {oo},  a  demand  function  6  :  V  — ►  R  such  that  J2veV  H^)  =  0>  ^^'^  ^  ^°^^  function 
c  :  £■  — »  R.  We  shall  use  the  notation  c{v,w)  —  —c{w,v)  for  v,w  EV  such  that  {w,v)  G  E.  In 
the  special  case  of  integer  demands,  we  shall  use  B  to  denote  the  maximum  absolute  value  of 
a  demand. 

A  pseudoflow  is  a  function  f  :  E  —^  R"*",  such  that  f{v,w)  <  u{v,w)  for  every  {v,w)  £  E. 
Given  a  pseudoflow  /,  we  define  the  excess  function  e/  :  V  — ►  R  by 

^ji'")=       H      f{w,v)-       ^      f{v,w)-b{v), 

w.(w,v)^E  w.(v,w)&E 

the  amount  by  which  the  net  flow  into  v  exceeds  the  demand.  We  say  that  a  node  v  has  excess 
if  ej{v)  is  positive,  and  has  deficit  if  it  is  negative.  For  a  subset  S  of  the  nodes  we  shall  use 
ef{S)  to  denote  X]ue5  ^/('^)-   ^'^^  ^  node  v,  we  define  the  ^ou;  conservation  constraint  by 


ef{v)  =  0     (flow  conservation  constraint).  (1) 

A  pseudoflow  /  is  a  transshipment  if  it  satisfies  the  flow  conservation  constraints  at  every 
node. 

The  residual  graph  with  respect  to  a  pseudoflow  /  is  given  by  Gj  =  (V,  £"/),  where  Ej  = 
{{v,w)\f{v,w)  <  u{v,w)  or  f{w,v)  >  0}. 

The  cost  of  a  pseudoflow  /  is  given  by 

c{f)  -      ^     f{v,w)c{v,w). 

The  transshipment  problem  is  that  of  finding  a  minimum-cost  [optimal)  transshipment  in  an 
input  network  (G,u,6,c).  In  the  uncapacttated  transshipment  problem  all  capacities  are  equal 

to  oo. 

In  order  to  simplify  the  presentation  we  restrict  our  attention  to  the  uncapacitated  trans- 
shipment problem  and  assume  that  the  graph  G  is  strongly  connected.  The  adaptation  of  the 
presented  results  to  the  general  transshipment  problem  is  straightforward.  Given  an  instance 
of  the  transshipment  problem  with  capacities,  we  can  construct  an  equivalent  uncapacitated  in- 
stance by  introducing  a  new  node  in  the  middle  of  every  arc.  The  dual  simplex  method  applied 
to  the  resulting  uncapacitated  problem  is  the  same  as  the  dual  network  simplex  algorithm  for 
the  original  problem.  This  construction,  however,  increases  the  number  of  nodes  to  0{m  +  n). 
The  assumption  that  the  graph  G  is  strongly  connected  can  be  satisfied  by  introducing  at  most 
2n  additional  expensive  arcs  {s,v)  and  {v,s)  for  some  node  s  and  every  v.  If  the  original  prob- 
lem is  feasible  then  no  optimal  solution  uses  the  additional  arcs.  This  assumption  guarantees 
that  the  transshipment  problem  is  feaisible  and  implies  the  following  characterization  of  the 
existence  of  an  optimal  solution. 

Theorem  2.1    There  exists  a  minimum-cost  transshipment  if  and  only  if  the  input  network 
contains  no  negative-cost  cycles. 

Linear  programming  duality  theory  provides  a  criterion  for  the  optimality  of  a  transship- 
ment. To  state  the  criterion  we  need  the  notions  of  a  price  function  and  a  reduced  cost  function. 
A  price  function  is  a  node  labelling  p  :  V  — ►  R.  The  reduced  cost  function  with  respect  to  a  price 
function  p  is  defined  by  Cp{v,  w)  —  c(ii,  w)  4-  p(i')  —  p('i')-  For  the  uncapacitated  transshipment 
problem  a  feasible  dual  solution  is  a  set  of  prices  such  that  each  arc  in  E  has  non-negative 
reduced  cost. 


Theorem  2.2   [5]  A  flow  f  is  an  optimal  solution  for  the  uncapacitated  transshipment  problem 
if  and  only  if  there  is  a  price  function  p  such  that,  for  each  arc  {v ,w)  G  E, 


procedure  Pivot(T,  (t;',  u;')); 

We  assume  that: 

T  is  a  dual  feasible  tree; 

p  is  the  price  function  that  corresponds  to  T; 

(v',w')  an  arc  such  that  f^(v',w')  <  0; 

Let  {w,  v)  E.  E  he  the  arc  with  minimum  reduced  cost  6  leaving  Htxv'.w'Yi 
Replace  the  arc  (v',w')  in  T  by  {v,w); 
Decrease  p  on  all  nodes  in  Hr,(v>.w')  by  <5; 
end. 


Figure  1:  A  Pivot  Step  of  the  Dual  Network  Simplex  Method. 

c.p{v,  w)  >Q  (dual  feasibility  constraints)      and  .  . 

Cp{y,  iz;)  >  0  =>  /(i;,  u;)  =  0     (complementarity  slackness  conditions). 

A  basis  of  the  linear  program  corresponding  to  the  uncapacitated  transshipment  problem 
is  the  set  of  columns  corresponding  to  the  arcs  in  a  spanning  tree  T  of  G.  A  transshipment  is 
a  basic  primal  solution  if  and  only  if  the  arcs  with  non-zero  flow  form  a  forest.  A  set  of  dual 
feasible  prices  is  a  basic  dual  solution  if  and  only  if  the  arcs  with  zero  reduced  cost  form  a 
spanning  subgraph. 

Let  T  be  a  spanning  tree  of  the  underlying  undirected  graph.  Consider  the  cut  obtained  by 
deleting  an  arc  {v,  w)  £  T  from  T.  Let  //r.fu.u/)  denote  the  set  of  nodes  that  are  on  the  same 
side  of  this  cut  as  w,  the  head  of  the  arc  (v,  w). 

For  a  tree  T  let  f^  denote  the  corresponding  bcisic  flow,  that  is, 

Similarly,  we  let  p^  denote  the  prices  defined  by  the  spanning  tree  T,  i.e.,  prices  such  that 
the  reduced  costs  of  the  arcs  in  the  tree  T  with  respect  to  prices  p     are  0. 

Let  p  =  p^  denote  the  prices  defined  by  the  tree  T.  The  spanning  tree  T  is  dual  feasible 
if  Cp{v,w)  >  0  for  all  arcs  {v,w)  E  E.  A  spanning  tree  T  is  called  primal  feasible  if  /  >  0. 
A  spanning  tree  is  optimal  if  it  is  both  prmial  feasible  and  dual  feasible.  Theorem  2.2  implies 
that  in  this  case  the  defined  flow  f^  is  an  optimal  transshipment. 

The  dual  network  simplex  algorithm  (see  Figure  1)  maintains  a  dual  feasible  tree  T.  An 
initial  dual  feasible  tree  can  be  found  by  a  shortest  path  computation.  In  the  case  that  all  arcs 
are  non-negative,  one  can  find  the  initial  tree  T*^  in  0{n)  pivots  [4];  in  case  the  costs  may  be 
negative,  one  can  find  T°  in  0{n^)  pivots  using  the  dual  simplex  algorithm  of  Balinski  [2].  A 
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pivot  step  of  the  dual  network  simplex  algorithm  can  be  applied  to  an  arc  {v' ,  w')  in  T  with 
f^{v' ,  w')  <  0.  The  pivot  step  changes  T  by  deleting  the  arc  (n',  w')  and  replacing  it  by  one 
of  the  arcs  {v,w)  of  minimum  reduced  cost  leaving  the  set  H-]-(v',w')-  The  arc  {v',w')  is  the 
leaving  arc,  and  {v,w)  is  the  entering  arc. 

Let  T'  =  T  —  {v' ,  w')  +  (u,  w)  be  the  tree  obtained  after  the  pivot.  It  is  easy  to  show  that  T' 
is  dual  feasible.  Moreover,  the  dual  objective  function  J2v^v  P(^)K^)  has  not  decreased  after 
the  pivot.  Let  6  be  the  reduced  cost  of  the  entering  arc  [v,  w).  We  may  obtain  the  price  vector 
p^  from  the  previous  price  vector  p^  by  subtracting  6  from  the  price  of  each  node  in  Hj(^,^^i\ 
and  keeping  all  other  prices  the  same. 

The  dual  simplex  algorithm  iterates  this  basic  pivot  step  until  the  current  tree  becomes 
primal  feasible.  Pivoting  strategies  give  rules  for  choosing  among  the  possible  leaving  and 
entering  arcs. 


3      A  simple  polynomial  time  dual  simplex  algorithm 

In  this  section  we  describe  how  to  use  an  excess  scaling  algorithm  to  guide  the  selection  of  dual 
pivots.  In  the  first  subsection  we  present  a  basic  subroutine  consisting  of  a  sequence  of  at  most 
n  pivot  steps  that  makes  the  next  augmentation  possible.  Next  we  show  how  to  use  a  version 
of  the  Edmonds-Karp  capacity  scaling  algorithm  to  derive  a  polynomial  dual  network  simplex 
algorithm  that  makes  at  most  O(n'^lognfi)  pivots.  Finally  we  give  a  strongly  polynomial 
version  with  an  0{n^  log  n)  bound  on  the  number  of  pivots.  The  algorithm  is  a  simplification  of 
the  polynomial  time  dual  simplex  algorithm  presented  in  [9].  Also,  as  opposed  to  the  algorithm 
presented  in  [9],  the  strongly  polynomial  algorithm  described  in  this  section  uses  only  the  usual 
arithmetic  operations  (in  fact,  only  additions  and  comparisons). 

We  shall  use  r  to  denote  the  special  root  node.  We  shall  use  T  to  denote  a  tree  rooted  at  r, 
and  use  Ty  to  denote  the  subtree  of  T  routed  at  node  v.  We  use  pred{v)  to  denote  u's  parent 
in  the  tree.  An  arc  (u,  w)  ^T  \s  called  downward  if  v  is  the  parent  of  w  in  the  tree  (  i.e.  {v,  w) 
points  away  from  the  root  node),  otherwise  (i;,  u;)  is  called  upward. 

The  algorithm  will  maintain  a  pair  (T,  /)  of  a  tree  and  a  pseudoflow  satisfying  the  following 
conditions: 

PI.  T  is  a  dual  feasible  tree. 

P2.  f(v,w)  >  Ofor  {v,w)  eT. 

P3.  f[v,w)  =  Ofor  {v,w)  ^T. 

P4.  ef{v)  >  0  for  each  v  ^  r. 

Recall  that  we  have  defined  /     as  the  bcisic  flow  corresponding  to  tree  T. 

Lemma  3.1  If{v,w)  is  a  downward  arc  of  T ,  then  f  {v,w)  =  f{v,w)  —  ef{Tu,).  If{v,w)  is 
an  upward  arc  of  T ,  then  f^{v,  w)  =  f[v,  w)  +  ej{Ty). 


procedure  MaKE-Good(T, /); 

We  assume  that: 

T  is  a  dual  feasible  tree; 

Let  p  be  the  price  function  that  corresponds  to  T; 

while  T  has  bad  arcs  do  begin 

let  S  denote  the  set  of  bad  nodes  for  T  w.r 

.t.  /; 

let  {v,w)  be  the  minimum  reduced  cost  arc 

leaving  5; 

let  6  *—  Cp(v,w); 

let  {v',  w')  be  the  first  bad  arc  on  the  path 

from  r  to  v; 

letT  ^T  +  (v,w)-{v',w')■ 

Decvease  p  on  all  nodes  in  Hx,(v',w')  by  6\ 

end; 

end. 

Figure  2:  Procedure  Make-Good. 

Proof:  One  can  obtain  the  flow  /^  from  the  flow  /  by  sending  ej{v)  units  of  flow  from  u  to  r 
for  each  node  v  ^  r.  The  increase  of  the  flow  in  an  upward  arc  {v,w)  is  ef{Tui)-  The  decrease 
of  the  flow  in  a  downward  arc  {v,w)  is  ef^T^).       | 

Suppose  that  the  pair  (/,  T)  satisfies  P1-P4.  An  arc  {v,  w)  in  T  is  called  badif  it  is  downward 
and  its  flow  is  0.  Otherwise,  it  is  called  good.  A  node  t;  G  T  is  called  good  if  every  arc  on  the 
path  from  r  to  u  is  good.  Otherwise,  node  v  is  called  bad. 

Corollary  3.2  Suppose  that  {f,T)  satisfies  Pl-P^-  Then  every  bad  arc  is  an  eligible  exiting 
arc. 


3.1     The  procedure  Make-Good 

We  will  now  describe  a  pivoting  procedure  Make-Good  (see  Figure  2)  for  transforming  a  tree 
T  satisfying  P1-P4  that  hzis  some  bad  arcs,  into  a  tree  in  which  all  arcs  are  good.  Given  a  tree 
T  and  flow  /,  Make-Good  proceeds  in  iterations.  In  each  iteration  it  considers  the  set  5  of 
bad  nodes  and  finds  a  minimum  reduced  cost  arc  (v,  w)  leaving  this  set  and  an  arc  {v',  w')  that 
is  the  first  bad  arc  on  the  path  from  r  to  v.  Then  it  removes  {v,w)  from  the  tree  and  adds 
{v',w')  instead.  Make-Good  terminates  when  there  are  no  more  bad  nodes. 

Since  {v',w')  is  a  bad  arc.  Corollary  32  implies  that  it  is  an  eligible  arc  for  exiting.  In 
general,  adding  the  arc  {v,w)  instead  of  {v',w')  does  not  maintain  dual  feasibility.  However, 
Make-Good  is  called  only  under  certain  conditions  that  ensure  that  it  executes  only  legal 
pivots. 

Make-Good  is  called  by  the  simplex  algorithms  described  in  the  subsequent  sections  after 
a  flow  augmentation  is  done  from  some  node  to  the  root  r,  if  this  augmentation  created  a  bad 

8 


Proof:  First  we  show  that  any  good  node  in  T  is  also  a  good  node  in  T' .  Suppose  that  t  is 
a  good  node  of  T.  Then  the  path  P  from  <  to  r  in  T  does  not  contain  the  bad  arc  {v',w')  and 
hence  P  is  also  a  path  in  T'.  It  follows  that  t  is  good  in  T'. 

New,  we  note  that  u;  is  a  good  node  in  T,  and  v  is  bad  in  T,  and  thus  v  becomes  good 
in  T'.  We  conclude  that  the  number  of  good  nodes  in  T'  is  greater  than  the  number  of  good 
nodes  in  T.       | 

Next  we  prove  that  this  sequence  of  pivots  can  be  performed  in  0{m  +  n  log  n)  time  using 
Fibonacci  heaps. 

Theorem  3.5  Suppose  that  tree  T  satisfies  PI-P4  and  all  of  the  bad  arcs  of  T  lie  on  some 
inverse-good  path  P  from  some  node  t  to  the  root  r.  Then  the  Make-Good  procedure  can  be 
implemented  in  0{m  +  nlogn)  time. 

Proof:  The  proof  is  based  on  a  closer  look  at  the  proof  of  Theorems  3.4  and  3.3.  The  algorithm 
will  maintain  the  first  bad  node  w'  on  the  path  from  r  tot,  and  for  every  bad  node  i;  it  maintains 
d{v)  =  minju^good)  Cp(i^,u')  and  the  good  node  Wy  on  which  the  minimum  is  attained. 

Recall  that  the  set  of  bad  nodes  is  the  set  of  nodes  in  T^i.  The  entering  edge  is  {v,  u;„)  for 
a  bad  node  v  with  d{v)  minimum.  The  leaving  edge  is  {pred{w'),w').  Next  we  have  to  show 
how  to  update  the  above  information.  The  new  w'  is  the  head  of  the  first  bad  arc  along  the 
path  from  v  to  t.  To  find  the  new  w'  after  the  update  takes  time  proportional  to  the  number 
of  nodes  from  v  to  the  new  w' .  Notice  that  these  are  the  nodes  that  become  good  during  this 
pivot. 

To  update  d{v)  for  the  remaining  bad  nodes  we  have  to  consider  all  edges  leaving  the  nodes 
that  became  good  during  this  pivot.  We  shall  maintain  d{v)  in  a  Fibonacci  heap.  Overall, 
there  are  at  most  m  updates  (one  for  every  arc).  The  total  time  spent  over  all  iterations  of 
Make-Good  on  searching  for  nodes  that  become  good  and  for  the  new  w'  is  bounded  by  0{n). 
The  claim  follows  since  updating  and  finding  the  minimum  can  be  done  in  0{m  +  n  log  n)  time. 
I 


3.2      A  polynomial  dual  simplex  algorithm. 

We  are  now  prepared  to  describe  the  polynomial  time  dual  simplex  algorithm.  The  algorithm 
finds  an  optimal  solution  for  the  minimum  cost  flow  problem  starting  with  the  shortest  path 
tree  T°  directed  from  node  r.  In  the  case  that  all  arcs  are  non-negative,  one  can  find  T  in 
0{n)  pivots  [4];  in  case  the  costs  may  be  negative,  one  can  find  T^  in  0{n^)  pivots  using  the 
dual  simplex  algorithm  of  Balinski  [2].  We  assume  that  T°  is  obtained  using  an  appropriate 
subroutine. 

In  our  algorithm,  we  will  keep  track  of  a  scaling  parameter  A  that  is  nondecreasing  from 
iteration  to  iteration.  We  will  refer  to  the  A-scaling  phase  as  all  iterations  in  which  the  param- 
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arc  by  reducing  to  zero  the  flow  on  one  of  the  downward  arcs.  Make-Good  transforms  the 
tree  into  a  new  tree  where  there  are  no  bad  arcs,  and  hence  we  can  do  an  augmentation  from 
any  node  with  positive  excess  to  the  root. 

Let  T  be  a  tree  and  let  /  denote  a  flow  satisfying  P1-P4.  We  call  a  path  P  in  T  from  some 
node  t  to  the  root  r  an  inverse-good  path  if  all  upward  arcs  on  P  have  positive  flow.  Notice 
that  if  we  were  to  reverse  the  direction  of  each  arc  on  an  inverse-good  path  then  the  resulting 
path  would  consist  of  good  arcs.  Observe  that  if  a  tree  had  no  bad  arcs,  an  augmentation  from 
some  node  to  the  root  can  introduce  bad  arcs  only  on  the  path  used  by  the  augmentation  and 
hence  all  bad  arcs  will  lie  on  an  inverse-good  path. 

Theorem  3.3  Suppose  that  the  pair  {T ,  f)  satisfies  PI-P4  and  all  of  the  bad  arcs  of  the  tree  T 
lie  on  some  inverse-good  path  P  from  some  node  t  to  the  root  r.  Then  the  steps  of  the  procedure 
Make-Good(T, /)  are  dual  simplex  pivots. 

Proof:  Let  5  denote  the  set  of  bad  nodes  of  T.  Let  {v,w)  denote  the  arc  that  is  pivoted  in  and 
let  (u',  w')  denote  the  arc  that  is  pivoted  out.  We  first  claim  that  S  is  the  set  of  descendants 
of  node  w' .  By  definition,  all  descendants  of  w'  are  bad.  To  see  the  converse,  suppose  that  s 
is  a  bad  node.  By  the  choice  of  {v',w'),  all  arcs  on  the  path  from  r  to  v'  are  good.  Let  P'  be 
the  path  from  s  to  r  and  let  {v" ,  w")  denote  a  bad  arc  of  P' .  By  hypothesis  {v" ,  w")  is  also  on 
path  P,  and  thus  s  is  a  descendent  of  w' . 

We  now  claim  that  the  first  pivot  of  procedure  Make-Good  is  a  dual  simplex  pivot.  By 
Corollary  3.2  the  arc  (t;',  w')  is  an  eligible  exiting  arc.  For  the  pivot  to  be  a  dual  simplex  pivot, 
the  entering  variable  must  be  the  least  cost  arc  from  the  subtree  T^i'  to  the  rest  of  the  graph. 
This  is  implied  by  the  choice  of  (w,  w),  and  the  fact  that  the  set  of  bad  nodes  of  T  is  the  set  of 
nodes  of  7^/. 

Let  T  be  the  tree  obtained  after  a  pivot  from  tree  T.  We  will  show  that  all  of  the  bad  arcs 
of  T'  lie  on  an  inverse  good  path  from  <  to  r  in  T'.  The  theorem  then  follows  by  induction. 

Node  i;  is  a  bad  node  of  T,  and  thus  the  path  from  v  to  r  intersects  the  path  P  at  some 
node,  say  node  s.  Since  all  bad  arcs  lie  on  P,  the  path  P'  from  w  to  s  is  good.  Let  P\  denote 
the  subpath  of  P  from  s  to  w\  and  let  Pj  denote  the  subpath  of  P  from  t  to  s.  Observe  that 
each  arc  of  T  \  Pj  \  P'  —  (t;',  u;')  has  the  same  direction  (upward  or  downward)  in  T'  as  it  does 
in  T.  Each  arc  of  Pi  and  P'  has  the  opposite  direction  in  T'  as  it  does  in  T.  It  follows  that 
Pi  is  transformed  by  the  pivot  from  an  inverse-good  path  into  a  good  path,  and  that  P  is 
transformed  from  a  good  path  to  an  inverse-good  path.  Subsequently,  all  bad  arcs  of  T'  lie  on 
P2  or  on  P',  and  thus  on  the  inverse  good  path  P2,  P'  from  <  to  r  in  T'.       | 

The  following  theorem  bounds  the  number  of  pivots  that  can  be  made  by  the  procedure 
Make-Good. 

Theorem  3.4  Suppose  that  the  pair  (T,  /)  satisfies  PI-P4  and  all  of  the  bad  arcs  of  the  tree  T 
lie  on  some  inverse-good  path  P  from  some  node  t  to  the  root  r.  The  procedure  Make-Good 
terminates  with  a  good  tree  T  in  at  most  n  —  1  pivots. 
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Proof:  First  we  show  that  any  good  node  in  T  is  also  a  good  node  in  T'.  Suppose  that  t  is 
a  good  node  of  T.  Then  the  path  P  from  <  to  r  in  T  does  not  contain  the  bad  arc  (v' ,  w')  and 
hence  P  is  also  a  path  in  T'.  It  follows  that  t  is  good  in  T' . 

Now,  we  note  that  if  is  a  good  node  in  T,  and  i;  is  bad  in  T,  and  thus  v  becomes  good 
in  T' .  We  conclude  that  the  number  of  good  nodes  in  T'  is  greater  than  the  number  of  good 
nodes  in  T.       | 

Next  we  prove  that  this  sequence  of  pivots  can  be  performed  in  0{m  +  n  logn)  time  using 
Fibonacci  heaps. 

Theorem  3.5  Suppose  that  tree  T  satisfies  PI-P4  and  all  of  the  bad  arcs  of  T  lie  on  some 
inverse-good  path  P  from  some  node  t  to  the  root  r.  Then  the  Make-Good  procedure  can  be 
implemented  in  0{m  +  rj  log  n)  time. 

Proof:  The  proof  is  based  on  a  closer  look  at  the  proof  of  Theorems  3.4  and  3.3.  The  algorithm 
will  maintain  the  first  bad  node  w'  on  the  path  from  r  to  t,  and  for  every  bad  node  v  it  maintains 
d{v)  =  min(^goo<j)  Cp{v,w)  and  the  good  node  Wy  on  which  the  minimum  is  attained. 

Recall  that  the  set  of  bad  nodes  is  the  set  of  nodes  in  T^i .  The  entering  edge  is  {v,  Wy)  for 
a  bad  node  v  with  d{v)  minimum.  The  leaving  edge  is  {pred{iu'),w').  Next  we  have  to  show 
how  to  update  the  above  information.  The  new  w'  is  the  head  of  the  first  bad  arc  along  the 
path  from  t;  to  t.  To  find  the  new  w'  after  the  update  takes  time  proportional  to  the  number 
of  nodes  from  v  to  the  new  w' .  Notice  that  these  are  the  nodes  that  become  good  during  this 
pivot. 

To  update  d{v)  for  the  remaining  bad  nodes  we  have  to  consider  ail  edges  leaving  the  nodes 
that  became  good  during  this  pivot.  We  shall  maintain  d{v)  in  a  Fibonacci  heap.  Overall, 
there  are  at  most  m  updates  (one  for  every  arc).  The  total  time  spent  over  all  iterations  of 
Make-Good  on  searching  for  nodes  that  become  good  and  for  the  new  w'  is  bounded  by  0{n). 
The  claim  follows  since  updating  and  finding  the  minimum  can  be  done  in  0{m  +  n  log  n)  time. 
I 


3.2     A  polynomial  dual  simplex  algorithm. 

We  are  now  prepared  to  describe  the  polynomial  time  dual  simplex  algorithm.  The  algorithm 
finds  an  optimal  solution  for  the  minimum  cost  flow  problem  starting  with  the  shortest  path 
tree  T  directed  from  node  r.  In  the  Ccise  that  all  arcs  are  non-negative,  one  can  find  T  in 
0{n)  pivots  [4];  in  case  the  costs  may  be  negative,  one  can  find  T^  in  0{n^)  pivots  using  the 
dual  simplex  algorithm  of  Balinski  [2].  We  assume  that  T°  is  obtained  using  an  appropriate 
subroutine. 

In  our  algorithm,  we  will  keep  track  of  a  scahng  parameter  A  that  is  nondecrejising  from 
iteration  to  iteration.  We  will  refer  to  the  A- scaling  phase  as  all  iterations  in  which  the  param- 
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procedure  ScAUNG-SlMPLEX(T); 

We  assume  that: 

T  is  a  dual  feasible  tree  and  all  arcs  are  good  subject  to  the  flow  /  =  0; 

Define  /  by  sending  A  units  of  flow  from  r  to  every  node  v  ^  r  in  T . 
while  A  >  l/(2n)  do  begin 
while  S/(A)  ^  0  do  begin 
let  n  be  a  node  in  5/ (A); 
send  A  units  of  flow  from  v  to  r  in  T. 
if  there  is  a  bad  arc  in  T 
then  begin 

call  Make-Good(T,/) 
end; 
end; 

A  *- A/2 
end; 
end. 


Figure  3:  Scaling-Simplex  algorithm. 

eter  A  has  a  fixed  value.  Throughout  the  algorithm,  the  tree  T  and  the  flow  /  will  satisfy  P5 
in  addition  to  P1-P4. 

P5.  All  flows  are  multiples  of  A. 

Let  5/(A)  =  {v  :  ef{v)  >  A}.  The  termination  of  the  algorithm  will  be  guaranteed  by  the 
following  lemma. 

Lemma  3.6   Suppose  that  the  supplies  and  demands  are  integral  and  T  is  dual  feasible.  Suppose 
that  f  is  a  nonnegative  flow  such  that 

•  f{v,  w)  =  0  for  every  {v,w)  ^  T, 

•  0  <  ef{v)  <  \/n  for  every  i;  ^  r. 

Then  T  is  primal  feasible. 

Proof:   By  Lemma  3.1,  f^{v,w)  >  f{v,w)  -  1.  The  flow  /^  is  integral,  and  /  is  non-negative. 
It  follows  that  /     is  non-negative.       | 

Procedure  Scaling-Simplex  is  described  in  Figure  3.  It  starts  with  an  initial  dual-feasible 
tree  T°  where  all  arcs  are  pointing  away  from  the  root,  and  a  parameter  A  =  2l'°8'    "•"  ".  The 
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initial  flow  is  constructed  by  sending  A  units  of  flow  from  the  root  to  each  one  of  the  nodes 
through  the  edges  of  the  tree.  The  procedure  proceeds  in  phases.  Each  phase  considers  the  set 
^/(A)  of  nodes  whose  excess  is  above  A,  and  iteratively  augments  a  flow  from  one  of  the  nodes 
in  this  set  to  the  root  by  A.  If  such  an  augmentation  creates  a  bad  arc,  the  tree  is  updated 
by  caUing  the  procedure  Make-Good  described  above.  A  phaise  ends  when  the  set  5/(A) 
becomes  empty.  Then  A  is  halved  and  a  new  phcise  is  started.  ScALiNG-SlMPLEX  terminates 
when  A  falls  below  l/(2n). 

Theorem  3.7  The  ScALlNG-SlMPLEX  algorithm,  if  started  from  a  shortest  path  tree  rooted 
at  r,  performs  dual  simplex  pivots,  and  finds  the  optimal  spanning  tree  for  the  uncapacitated 
transshipment  problem  after  0{n'^  log  nB)  pivots,  and  in  0(n(m  +  n  log  n)  log  nfl)  time. 

Proof:  We  claim  that  the  algorithm  maintains  properties  P1-P5.  First  note  that  all  flows  are 
multiples  of  A  throughout  the  algorithm,  and  so  P5  is  satisfied.  We  now  claim  that  all  excesses 
are  strictly  positive.  It  is  true  initially,  since  initially  A  >  B  by  definition.  During  the  algorithm 
we  send  A  units  of  flow  only  from  nodes  with  excess  more  than  A,  and  so  all  excesses  remain 
positive  after  an  augmentation.  Therefore  P4  is  satisfied. 

Since  all  arc  flows  are  multiples  of  A,  one  can  send  A  units  of  flow  on  any  good  path.  For 
this  reason,  P2  and  P3  are  satisfied  throughout.  Initially  T  is  a  dual  feasible  tree  with  all 
arcs  being  good.  After  sending  A  units  of  flow  on  a  good  path  P  from  some  node  v  to  r  the 
path  P  becomes  an  inverse-good  path.  If  P  has  any  bad  arcs,  then  the  procedure  Make-Good 
takes  0{n)  pivots  and  0{m  +  nlogn)  time  by  Theorems  3.4  and  3.5.  Property  PI  is  satisfied 
throughout  the  execution  of  the  procedure  by  Theorem  3.3. 

ScALlNG-SlMPLEX  terminates  with  A  <  l/(2n).  At  the  end  of  the  last  scaling  phase  we 
obtain  a  tree  T  and  a  flow  /.  At  this  point,  ef(v)  <  \/n  for  each  node  w  /  r.  Lemma  3.6 
implies  that  T  is  primal  feasible,  and  hence  T  is  an  optimal  basis. 

During  a  scaling  phase  each  node  in  n  ^  r  starts  with  an  excess  of  less  than  2A  and  ends 
with  a  positive  excess.  Each  augmentation  reduces  the  excess  of  a  node  d  ^  r  by  A.  Thus 
the  number  of  augmentations  per  scaling  phase  is  at  most  n.  The  algorithm  terminates  when 
A  <  l/(2n),  therefore  the  total  number  of  augmentations  is  0{n\ognB). 

The  time  between  successive  augmentations  is  0{m  -\-  nlogn)  by  Theorem  3.5,  and  the 
number  of  pivots  is  0(n)  by  Theorem  3.4.  We  can  conclude  that  the  total  number  of  pivots  is 
0{n^  log  nB)  and  the  total  running  time  is  0{n[m  +  n  log  n)  log  nB).       | 


3.3      A  strongly  polynomial  dual  simplex  algorithm. 

We  will  modify  the  dual  simplex  algorithm  described  in  the  previous  section  to  make  it  strongly 
polynomial.  We  will  divide  scaling  phases  into  two  types,  according  to  whether  a  pivot  was 
performed  during  the  phase,  and  will  bound  the  number  of  scaling  phases  that  involve  pivots 
by  0{n  logn).  The  second  type  of  scaling  phase  does  not  involve  pivots  and  hence  the  spanning 
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tree  at  the  beginning  of  the  scahng  phase  is  the  same  as  the  spanning  tree  at  the  end  of  the 
scaling  phase.  The  number  of  such  scahng  phases  for  the  Scaling-Simplex  algorithm  is  at 
least  lognfl.  We  will  show  how  to  modify  the  algorithm  appropriately  so  that  the  number  of 
these  scaling  phases  is  0{n). 

We  first  show  that  the  number  of  scaling  phases  in  which  some  pivot  takes  place  is  0{n  log  n). 
To  prove  this  bound,  we  first  observe  that  all  flow  changes  in  the  algorithm  are  due  to  sending 
flow  on  a  path  to  the  root.  Since  the  root  node  is  the  only  node  with  a  negative  excess 
throughout  the  algorithm,  the  total  flow  change  on  any  arc  is  bounded  at  —ef{r).  Thus,  we 
have  the  following  lemma. 

Lemma  3.8  Suppose  that  ef{r)  <  0  throughout  the  algorithm,  and  that  each  flow  change  is 
the  result  of  sending  flow  from  some  node  v  :^  r  to  node  r.  If  f{v,w)  >  — e/(r)  for  some  flow 
f  obtained  in  the  algorithm,  then  f'{v,w)  >  — ey(r)  for  all  subsequent  flows  f  obtained  by  the 
algorithm,  including  the  optimum  flow. 

We  will  refer  to  an  arc  {v,w)  as  strongly  feasible  if  f{v,w)  >  — ey(r).  The  above  lemma 
implies  that  once  an  arc  becomes  strongly  feasible  it  will  stay  strongly  feasible  throughout  the 
rest  of  the  algorithm,  and  it  will  have  a  strictly  positive  flow  in  the  optimum  solution  obtained 
by  the  algorithm. 

The  reason  that  we  wanted  arc  flows  to  be  multiples  of  A  was  so  that  we  could  be  assured 
that  we  could  send  A  units  of  flow  on  any  good  path.  The  above  lemma  implies  that  it  is 
sufficient  to  weaken  property  Pb  and  require  that  the  flow  will  be  multiple  of  A  on  arcs  that 
are  not  strongly  feasible.  We  will  also  need  to  strengthen  property  P4. 

P'4.   A/2  <  ej{v)  <  3A  for  each  v^r. 

P'5.   For  each  {v,uj)  for  which  f{v,w)  <  -ey(r),  f{v,w)  is  an  integral  multiple  of  A. 

For  convenience,  we  will  write  P'l-P'5  to  mean  PI,  P2,  P3,  P'4,  P'5. 

Lemma  3.9  //(T,/,  A)  satisfies  P'l-P'5,  and  if  v  is  a  good  node,  then  one  can  send  A  units 
of  flow  from  v  to  r,  and  the  resulting  flow  f  is  non-negative. 

The  strengthening  of  property  P4  is  needed  for  the  next  lemma. 

Lemma  3.10  Suppose  that  the  triple  (T,/,  A)  satisfies  P'l-P'5.  Suppose  further  that  {v,w)  is 
a  bad  arc  for  spanning  tree  T.  Then  within  an  additional  Z  +  '2\\ogn\  scaling  phases  the  number 
of  strongly  feasible  arcs  will  increase. 

Proof:  Since  arc  {v,w)  is  bad,  it  is  a  downward  arc  and  f{v,w)  =  0.  Therefore,  /  is  zero  on 
all  arcs  leaving  T^,.   By  property  P'4,  each  node  in  T^  has  excess  at  least  A/2.   Now  let  A'  be 
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the  scaling  factor  3  +  2nogn]  scaling  phaises  later,  and  let  /'  denote  a  flow  after  the  A'  scaling 
phase.  Then  A'  <  A/(8n2).  By  property  P'4,  ej{T^)  >  \T^\A/2  >  4n2|T^|A'.  Moreover, 
ef'{Tw)  <  3\Tu,\A'.  Thus,  at  least  3n2|Tu;|A'  units  of  flow  have  been  sent  from  T^  to  V  \  T^, 
while  transforming  flow  /  into  flow  /'.  At  most  n  -  1  arcs  have  positive  flow  in  /'.  This  implies 
that  there  must  be  some  arc  leaving  Tw  with  flow  greater  than  3n|Tu;|A'  >  3nA'  >  -e//(r), 
where  the  last  inequality  follows  from  P'4.  By  definition,  such  an  arc  is  strongly  feasible.  Since 
there  are  no  arcs  leaving  T^  with  positive  flow  subject  to  /,  it  follows  that  the  number  of 
strongly  feasible  arcs  has  increased.       | 

We  have  shown  that  if  some  arc  becomes  bad  during  a  scaling  phase,  then  within  3  +  2|'log  n] 
additional  scaling  phases  there  is  a  new  strongly  feasible  arc.  Since  there  are  at  most  n  arcs 
that  can  become  strongly  feasible,  it  follows  that  there  are  0(n  log  n)  iterations  in  which  an  arc 
becomes  bad,  causing  the  algorithm  to  execute  at  lecist  one  pivot. 

Corollary  3.11    The  number  of  scaling  phases  in  which  some  pivot  takes  place  is  O(nlogn). 

The  above  corollary  implies  a  strongly  polynomial  bound  on  the  number  of  pivots.  In  order 
to  show  a  strongly  polynomial  bound  on  the  running  time  we  have  to  limit  the  number  of 
scaling  phases  that  do  not  execute  pivots.  Let  T  be  the  current  spanning  tree  and  let  /  denote 
the  current  flow  at  the  A  scaling  phase.  Suppose  that  (T,/,  A)  satisfies  P'l  -  P'5.  We  will 
determine  a  new  scaling  factor  A'  <  A  and  a  flow  /',  such  that  (T,/',  A')  also  satisfies  P'l  - 
P'5,  and  such  that  there  is  no  flow  /"  for  which  (r,/",A/8)  satisfies  P'l-P'5.  Consequently, 
if  we  continue  the  scaling  algorithm  starting  with  flow  /',  then  within  3  scaling  phases  the 
spanning  tree  T  cannot  be  feasible,  and  there  has  to  be  a  pivot. 

Suppose  that  (T,/,A)  satisfies  P'l  -  P'5.  Let  A'  =  max(-2/^(u,  u;)/|T(i<;)|  :  {v,w)  €  T). 
Observe  that  A'  <  0  if  and  only  if  T  is  an  optimal  tree.  We  will  show  subsequently  that  there 
is  a  flow  /'  such  that  (T,/',  A')  satisfies  P'l-P'5.  First  we  show  that  there  can  be  no  flow  /" 
such  that  (T,/",  A'/S)  satisfies  P'l-P'5. 

Lemma  3.12   //A'  >  0  then  there  is  no  flow  f"  such  that  (T,/",  A'/S)  satisfies  P'l-P'5. 

Proof:  Let  /"  be  such  a  flow.  Consider  the  flow  on  arc  {v,  w)  where  A'  =  —2f^{v,  w)/\T{w)\. 
By  Lemma  3.1  and  Property  P'4,  f"{v,  w)  =  f'^{v,  w)  +  ejn{T^)  =  -A'|r(tx;)|/2  +  SfiT^)  < 
-A'\T{w)\/2  +  |T(k;)|(3A78)  <  0.       I 

Next  we  show  how  to  find  a  feasible  flow  /'  such  that  (T,  /',  A')  satisfies  P'l-P'5,  assuming 
that  A'  <  A/3.  The  procedure  Make-Flow  is  given  in  Figure  4.  We  start  by  creating  an 
excess  of  3A'/2  at  each  node  in  t;  /  r.  Let  g  denote  the  resulting  flow.  We  then  examine  nodes 
in  the  reverse  of  a  breadth  first  search  ordering  (i.e.,  we  start  at  the  leaves  and  work  towards 
the  root),  and  for  each  node  v  examined  we  modify  g  to  make  sure  that  the  flow  on  {pred{v),  v) 
satisfies  P'5. 
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procedure  MaKE-Flow(T); 

A'  -  maX(„,^)eT(-2/^(t^,  w)/\T(w)\). 

Define  g  by  sending  6(t;)  +  3A'/2  units  of  flow  from  r  to  every  node  v  ^  r  in  T. 
for  i»  G  V  in  a  the  reverse  BFS  ordering  do  begin 
Let  w  =  pred(v). 

if  (w,v)  is  a  downward  arc  and  (w,v)  is  not  strongly  feasible 
then  begin 

send  f(w,  v)  mod  A'  units  of  flow  from  v  to  node  r  in  T. 
end; 
if  {v,w)  is  an  upward  arc  and  {v,w)  is  not  strongly  feasible 
then  begin 

send  —f{w,v)modA'  units  of  flow  from  v  to  node  r  in  T. 
end; 
end; 
end. 


Figure  4:  Make-Flow  algorithm. 

Lemma  3.13  Suppose  that  A'  <  A/3.  Then  the  flow  g  computed  by  Make-Flow  satisfies 
eg{v)  <  ej{v)  for  every  v  ^  r. 

Proof:  By  property  P'4,  we  have  that  ej{v)  >  A/2,  A'  <  A/3  and  eg{v)  <  3A'/2.  This  implies 
the  lemma.       | 

Lemma  3.14  Suppose  that  (T,/,  A)  satisfies  P'l  -  P'5,  and  A'  <  A/3.  Then  {T,g,A')  com- 
puted by  Make-Flow  satisfies  P'l-P'5  and  all  arcs  of  T  are  good. 

Proof:  PI  and  P3  are  satisfied  by  definition.  By  construction,  g  satisfies  P'4  and  P'5.  It 
remains  to  show  that  the  resulting  flow  g{v,uj)  is  non-negative  on  all  upward  arcs  in  T  and 
positive  on  all  downward  arcs.  This  will  imply  imply  P2,  and  the  fact  that  all  arcs  are  good. 

Since  /  is  nonnegative  and,  by  Lemma  3.13,  we  have  eg{v)  <  ef{v)  for  v  ^  r,  the  flow  g 
on  upward  arcs  is  nonnegative.  Suppose  now  that  {v,xv)  is  a  downward  arc.  Then  g{v,w)  = 
f'^{v,w)  +  eg{T^)  >  -\T{w)\A.'/2  +  eg{J^)  >  0,  since  each  node  in  T^  has  at  least  A72  excess. 
I 

Notice  that  Lemma  3.13  implies  that  the  procedure  Make-Flow  can  be  also  thought  of  as 
modifying  /  by  sending  flow  to  the  root  r  from  some  other  nodes. 

Lemma  3.15  Suppose  that  (T,/,  A)  satisfies  P'l  -  P'5,  and  A'  <  A/3.  Then  every  arc  that 
2s  strongly  feasible  subject  to  f  is  also  strongly  feasible  subject  to  g. 

Proof:  The  lemma  is  implied  by  Lemmas  3. 13  and  3.8.       I 
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procedure  Strong-Scaling-Simplex(T); 

We  assume  that; 

T  is  a  dual  feasible  tree  and  all  arcs  are  good  subject  to  the  flow  /  =  0; 

A  <—  GO. 

while  T  is  not  optimum  do  begin 

let  A'  =  max(-2/^(t;,  w)/\T{w)\  :  (v,  w)  G  T) 
if  A'  =0 
then  begin 

quit  and  return  the  optimal  tree  T; 
end; 
if  A'  >  A/6 

then  A  ^  A/2 
else  begin 

call  Make-Flow(T),  /  «—  ff  and  A  «—  A'. 
end; 
while  5/ (3 A/2)  ^  0  do  begin 

let  V  be  a  node  in  5/(3A/2);  send  A  units  of  flow  from  v  to  r. 
if  there  is  a  bad  arc  in  T 
then  begin 

call  Make-Good(T,/) 
end; 
end; 
end; 
end. 


Figure  5:  SXRONG-ScALING-SlMPLEX  algorithm. 

Theorem  3.16  The  SxRONG-ScALING-SlMPLEX  algorithm,  if  started  from  a  shortest  path  tree 
rooted  at  r,  performs  dual  simplex  pivots,  and  finds  the  optimal  spanning  tree  for  the  uncapac- 
itated  transshipment  problem  after  0(n^  log  n)  pivots,  and  0{n^{m  +  nlogn)logn)  time. 

Proof:  The  tree  T  is  good  before  each  augmentation,  and  all  bad  arcs  are  on  an  inverse  good 
path  before  each  call  to  Make-Good.  The  fact  that  the  algorithm  maintains  P'l-P'5  can  be 
proved  analogously  to  the  proof  of  Theorem  3.7.  The  procedure  Make-Flow  takes  0{n)  steps. 
Subsequently,  within  0(1)  scaling  phcises  there  will  be  a  bad  arc,  and  by  Lemma  3.10  in  C>(log  n) 
scaling  phases  a  new  arc  will  become  strongly  feasible.  There  can  be  at  most  n  strongly  feasible 
arcs  at  a  time.  Therefore,  Lemma  3.15  implies  that  there  will  be  at  most  O(nlogTi)  scaling 
phases.  Now  Theorems  3.4  and  3.5  imply  the  bounds  on  the  number  of  pivots  and  the  running 
time.       I 

Using  the  reduction  from  the  minimum-cost  flow  problem  to  the  uncapacitated  transship- 
ment problem  we  obtain  the  following  corollary. 

Corollary  3.17    The  (capacitated)  transshipment  problem  can  be  solved  by  the  dual  network 
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simplex  method  in  0{Tn   log  n)  pivots,  and  the  algorithm  can  be  implemented  to  run  in  0{m'^{Tn+ 
n  log  n)  log  n)  time. 


4      An  Efficient  Pivoting  Strategy 

In  this  section  we  show  how  to  use  a  variation  on  Orlin's  [11]  fast  minimum-cost  transshipment 
algorithm  to  design  a  pivoting  strategy  that  decreases  the  number  of  pivots  for  the  dual  network 
simplex  algorithm  to  0(mn  log  n). 

The  simplex  algorithm  in  the  previous  section  can  be  viewed  as  an  implementation  of  the 
corresponding  scaling  algorithm.  During  the  scaling  algorithm  all  nodes,  except  the  root  r, 
have  positive  excess.  The  simplex  implementation  maintains  a  dual  feasible  tree  T  such  that 
an  augmentation  from  every  node  v  to  the  root  r  is  possible  in  the  tree.  Whenever  the  scaling 
algorithm  augments  the  flow  from  some  node  t;  to  r,  this  augmentation  can  be  done  in  the  tree. 

Such  a  direct  simplex  implementation  of  the  faster  scaling  algorithm  is  not  possible.  For 
some  pairs  of  nodes  that  the  scaling  algorithm  might  choose  for  an  augmentation  there  might  be 
no  sequence  of  pivot  steps  that  makes  an  augmentation  between  these  two  nodes  possible  in  the 
tree.  The  version  of  the  transshipment  algorithm  used  here  allows  greater  freedom  in  choosing 
the  augmentation.  We  show  that  a  simple  pivoting  strategy  can  be  viewed  as  implementing 
one  of  the  possible  choices.  In  eff"ect,  we  let  the  simplex  algorithm  choose  among  the  possible 
augmentations. 

In  Section  4.1  we  give  the  modified  version  of  Orlin's  scaling  algorithm,  and  review  the 
proof  of  its  running  time.  In  Section  4.2  we  use  this  algorithm  to  design  a  dual  network  simplex 
algorithm  that  takes  at  most  0(mn  log  n)  pivots. 

4.1     Modified  version  of  Orlin's  excess  scaling  algorithm 

The  main  idea  of  all  excess  scahng  algorithms  is  to  maintain  a  feasible  pseudoflow  and  a  price 
function,  such  that  the  reduced  cost  of  every  arc  is  non-negative,  and  the  reduced  cost  of  arcs 
with  positive  flow  is  zero.  The  initial  flow  is  zero,  initial  prices  can  be  computed  by  a  shortest 
path  computation.  We  repeatedly  augment  the  flow  along  shortest  paths  from  a  node  with 
positive  excess  to  a  node  with  negative  excess,  gradually  reducing  all  the  excesses  to  zero.  Two 
observations  justify  this  method:  (1)  moving  the  flow  along  a  minimum-cost  path  preserves  the 
invariant  that  the  current  pseudoflow  has  the  minimum-cost  among  all  the  pseudoflows  with 
the  same  excess;  (2)  a  shortest  path  computation  suffices  both  to  find  a  path  along  which  to 
augment  the  flow  and  to  find  appropriate  price  changes  that  preserve  the  nonnegativity  of  the 
reduced  costs.  Let  Augment(5,  T,  A)  denote  the  subroutine  that  augments  the  flow  by  A 
along  a  shortest  path  from  some  node  in  S  to  some  node  in  T. 

The  scaling  algorithm  that  was  the  basis  of  the  simplex  algorithm  in  the  previous  section 
maintains  that  all  nodes  except  the  root  r  has  positive  excess.  In  Orlin's  scaling  algorithm  we 
need  to  augment  the  flow  by  more  then  the  excess  at  the  end  on  the  path,  thereby  changing 
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the  sign  of  the  excess  during  the  algorithm.  Let  us  define  the  set  of  nodes  with  large  deficit 
Tf{A)  -  {v  £  V  :  ef{v)  <  -A}  analogously  to  5/(A). 

The  algorithm  maintains  a  scaling  parameter  A,  and  it  consists  of  a  number  of  scaling 
phases.  Each  phase  consists  of  a  sequence  of  augmentations.  The  main  difference  between  this 
algorithm  and  the  one  used  in  the  previous  section  is  that  at  the  A-scaling  phase  we  either 
augment  the  flow  by  A  from  some  node  in  5  =  5/(^^A)  to  a  node  in  T  =  Tf{^A),  or, 
reversing  the  roles  of  5  and  T,  from  some  node  in  5  =  5/(^A)  to  a  node  in  T  =  T/(^^A) 
until  both  5/(^^^^A)  and  Tj(^^A)  become  empty.  At  this  point  we  divide  A  by  two  and  a 
new  scaling  phase  starts. 

The  idea  of  making  this  algorithm  strongly  polynomial  is  similar  to  the  previous  one.  Notice 
that  the  sets  5/(^A)  and  r/(^A)  are  monotone  decreasing  during  the  A  scaling  phase. 
Therefore,  there  can  be  at  most  n  augmentations  per  scaling  phcise.  The  amount  of  flow  moved 
during  the  A  scaling  phcise  is  no  more  than  nA,  and  if  an  arc  has  more  than  2nA  flow  then  it 
will  have  positive  flow  in  all  subsequent  phases.  Such  arcs  will  be  declared  strongly  feasible. 

More  precisely,  the  algorithm  Fast-Excess-Scaling  maintains  a  set  of  strongly  feasible 
arcs  Es  and  a  pseudoflow  /  that  satisfies  the  following  properties. 

51.  /(w,  u;)  >  0  on  every  {v,  w)  ^  Eg. 

52.  f{v,  u;)  >  0  on  every  {v,w)  6  E,. 

53.  Cp{v,w)  >  0  for  aU  arcs  in  the  residual  graph. 

54.  At  most  one  node  has  non-zero  excess  in  every  connected  component  of  Gj  =  {V,  Eg). 

55.  f{v,w)  is  a  multiple  of  A  on  every  arc  {v,w)  ^  E^. 

An  arc  becomes  strongly  feasible  if  it  carries  at  least  5nA  flow.  The  algorithm  Fast- 
Excess-Scaling  is  shown  in  Figure  6.  Initially,  A  =  max^gv  1^(^)1-  Each  time  we  declare 
an  arc  strongly  feasible  we  do  a  special  augmentation  to  collect  the  excess  of  each  connected 
component  of  E,  to  a  single  node.  After  a  scaling  phase  A  is  divided  by  2  as  long  as  the  current 
pseudoflow  is  non-zero  on  some  arc  not  in  E,.  When  every  arc  with  non-zero  flow  is  in  Es  the 
scaling  is  restarted  by  setting  A  to  be  the  maximum  absolute  value  of  a  current  excess. 

Lemma  4.1   Algorithm  Fast-Excess-Scaling  maintains  properties  S1-S5. 

Proof:  Notice  that  if  S/(^A)  is  not  empty  then  T/(^A)  must  not  be  empty  since  excesses 
sum  to  zero.  Similarly  if  T/(^A)  is  not  empty  then  5/(^A)  must  not  be  empty.  This  proves 
that  the  augmentations  are  possible. 

Throughout  the  algorithm  we  shall  maintain  a  pseudoflow  /  and  a  price  function  p.  It  is 
easy  to  see  that  SI,  and  S3-S5  is  satisfied  throughout.  Strongly  feasible  arcs  have  reduced  cost  0 
by  properties  S2  and  S3,  therefore  augmentation  through  a  connected  component  of  G,  can  be 
done  through  the  strongly  feasible  arcs.  Note,  that  this  version  of  the  algorithm  can  move  more 
than  nA  flow  during  the  A-scaling  phase.  Each  time  an  arc  is  declared  strongly  feasible,  we  do 
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procedure  Fast-Excess-Scaling; 
A  «—  majCt,gv  \b{v)\; 
while  A  ^  0  do  begin 

while  SfC^A)  U  Tf(^A)  ^  Id  do  begin 
if5/(^A)_#0 
then  begin 

5^5;(^A); 
T^Tji^A); 
CaU  AuGMENT(5,r,A); 

Add  {v,  w)  to  E,  for  every  {v,  w)  such  that  /(v,  w)  >  5nA; 

Use  ^j  arcs  to  collect  the  excess  of  the  new  connected  component  of  G,  to  one  node; 
end; 

else  begin 
S^Sfi^A); 
T^Tfi^A); 
Call  AUGMENT(5,r,A); 

Add  (v,w)  to  Eg  for  every  {v,w)  such  that  f{v,w)  >  5nA; 

Use  E,  arcs  to  collect  the  excess  of  the  new  connected  component  of  G,  to  one  node; 
end; 
end; 
end; 

if  /  is  zero  on  all  arcs  not  in  E, 
then  A  «—  maxtigv  |6(f)|; 
else  A  =  A/2; 
end. 


Figure  6:  Algorithm  Fast-Excess-Scaling. 


a  special  augmentation  to  collect  the  excess  of  a  connected  component  of  Gs  This  might  yield 
another  nA  amount  of  flow  to  be  moved  through  strongly  feasible  arcs.  Collecting  the  excess  in 
a  connected  component  of  Ga  can  yield  an  increcised  excess  at  some  nodes,  and  therefore  extra 
augmentations.  The  number  of  augmentations  during  a  phase  is  at  most  n  plus  the  number  of 
contractions.  Therefore  the  overall  amount  of  flow  moved  during  and  after  the  A  scaling-phase 
is  at  most  5nA  and  arcs  that  carry  more  than  5nA  flow  will  never  become  empty,  and  hence 
S2  is  satisfied.       | 

It  will  help  the  analysis  to  consider  the  arcs  in  E3  as  contracted,  and  hence  the  connected 
components  of  {V,  £',)  are  pseudo-nodes  of  the  contracted  graph.  We  will  use  V  to  denote  node 
set  of  the  contracted  graph. 

The  main  idea  in  the  analysis  is  to  show  after  a  node  v  (either  in  V  or  a  pseudo-node) 
participated  as  an  end-point  of  a  shortest  path  computation,  after  at  most  O(logn)  additional 
scaling  phcises,  an  arc  adjacent  to  v  will  become  strongly  feasible.  For  a  node  d  in  V  this  follows 
from  the  fact  that  v  ^  5/(A)  UTf{A)  unless  |6(t;)|  >  A,  and  therefore  O(logn)  scaling  phases 
after  the  node  first  served  as  a  starting  node,  an  arc  incident  to  v  will  carry  enough  flow  to  be 
contracted.  However,  |6(u)|  >  |e^(t;)|  might  not  hold  for  pseudo-nodes. 
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We  say  that  a  node  v  is  active  during  a  scaling  phase  if  u  £  5/(^^A)  UT/(^^A)  at  some 
time  during  the  phase. 

Lemma  4.2  The  number  of  shortest  path  computations  during  a  phase  is  bounded  by  the  num- 
ber of  active  nodes  during  the  phase. 

Theorem  4.3   A  node  v  can  be  active  in  at  most  0(log  n)  phases  before  it  is  contracted. 

Proof:  A  pseudo-node  can  become  active  once  due  to  contraction.  However,  when  a  node  v  is 
active  for  the  second  time,  it  must  already  exist  at  the  end  of  the  previous  scaling  phase.  Let 
A  denote  the  scaling  parameter  in  the  phcise  when  v  is  active  for  the  second  time.  If  A  for  this 
phase  was  defined  by  maLX^^y  ef{w)  then  ef{v)  —  b{v).  Otherwise,  the  scaling  parameter  in 
the  previous  phase  is  2A.  At  the  end  of  the  2A-scaling  phase  both  5/(^2A)  and  r/(^2A) 
are  empty.  Therefore,  at  the  beginning  of  this  phase,  we  have  ^^A  <  |e/(t;)|  <  '"~  ^A. 
But  b{v)  -  ef{v)  is  an  integer  multiple  of  2A.  This  implies  that  |6(t;)|  >  :^A  >  ^e/(t;).  In 
either  case,  after  at  most  O(logn)  more  scahng  phases  the  scaling  parameter  A  will  be  less  than 
\b{v)\/{5n^).  At  the  end  of  that  scaling  phase  the  flow  /  will  satisfy  the  following  inequality. 

I  J2  /(«'  "^)l  =  IK^)  -  ^f{y)\  >  Hv)\  -  A  >  (5r»2  -  1)A.  (4) 

wev 

Consequently,  at  least  one  arc  incident  to  v  carries  more  than  5nA  flow,  and  hence  v  will  be 
contracted.       | 

Lemma  4.2  and  Theorem  4.3  bound  the  number  of  shortest  path  computations  during  the 
algorithm.  All  other  work  takes  linear  time  per  scaling  phase.  At  least  one  arc  is  contracted 
in  each  group  of  O(logn)  scaling  phases,  and  therefore,  there  are  at  most  O(nlogn)  scaling 
phases. 

Theorem  4.4  Algorithm  Fast-Excess-Scaling  solves  the  transshipment  problem  in  0{n  log  min{n,  B}) 
computations  of  single-source  shortest  paths  in  networks  with  non-negative  lengths. 

4.2      Decreasing  the  number  of  pivot  steps 

The  sequence  of  pivot  steps  in  the  dual  network  simplex  algorithm  in  the  previous  section  is 
guided  by  a  capacity  scaling  algorithm.  The  main  needed  change  in  the  scaling  algorithm  to 
guide  the  pivot  selection  is  that  all  augmentations  are  carried  out  in  the  tree.  In  this  section 
we  give  a  similar  simplex  implementation  of  algorithm  Fast-Excess-Scaling.  If  the  set  T  is 
not  reachable  from  S  over  residual  arcs  in  the  tree,  then  the  algorithm  initiates  a  sequence  of 
dual  simplex  pivot  steps,  changing  the  tree  T  into  one  in  which  an  augmentation  in  the  tree 
is  possible.    The  first  call  to  Augment  is  replaced  by  Tree-Augment+,  and  the  second  call 
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procedure  Tree-Augment+(5,T,  A,T,/); 

We  assume    that: 

T  is  a  dual  feasible  tree; 

p  is  the  price  function  which  corresponds  to  T; 

/  is  a  pseudoflow  that  is  zero  outside  of  the  tree  T. 

Choose  node  s  €  5; 

R  «—  Nodes  reachable  from  s  over  residual  arcs  in  T  subject  to  the  flow  /; 

while  R(T)  PI  T  =  0  and  there  is  an  arc  {v,  w)  entering  R{T)  such  that  f^ {v,  w)  <  0  do  begin 

{v,w)  ^—  an  arc  entering  R(T)  such  that  f~{v,w)  <  0; 

Call  PlVOT(T,(i;,u;)); 

Update  R(T)\ 
end; 

t  ^  A  node  in  R(T)nT; 

P  *—  A  path  from  s  to  <  in  through  residual  arcs  in  T; 
Move  A  units  of  flow  from  s  to  t  edong  P; 
end. 


Figure  7:  Pivot  Steps  that  Make  an  Augmentation  from  s  Possible. 

by  Tree- Augment".  The  Tree- Augment"*"  procedure  first  chooses  a  node  s  £  S,  and  then 
executes  a  sequence  of  pivot  steps  until  there  is  some  t  E  T  reachable  from  s  in  the  tree.  The 
Tree- Augment"  procedure  first  chooses  a  node  t  E  T,  and  then  execute  pivot  steps  until  there 
is  some  s  £  S  such  that  t  is  reachable  from  it  through  the  residual  arcs  of  the  tree.  Below 
we  show  that  the  properties  of  the  Fast-Excess-Scaling  algorithm  allow  us  to  prove  that  the 
number  of  pivots  needed  to  be  done  between  any  two  augmentations  is  small. 

The  procedure  Tree-Augment"*"  (see  Figure  7)  starts  by  choosing  a  node  s  E  S.  Let  R(T) 
denote  the  set  of  nodes  reachable  from  s  over  residual  arcs  in  the  tree  T.  If  the  set  R(T)  H  T 
is  not  empty,  then  we  augment  the  flow  along  the  path  in  the  tree  from  s  to  some  node  t  in 
the  intersection.  If  R(T)  fl  T  is  empty,  then  we  iteratively  choose  an  arc  {v,w)  entering  R{T) 
such  that  /■^(w,  w)  <  0,  and  pivot  on  this  arc,  until  the  intersection  R{T)  D  T  is  not  empty. 
Then  we  augment  the  flow  from  s  to  some  node  t  in  the  intersection.  Since  Tree-Augment" 
is  analogous,  we  omit  its  description. 

Next  we  need  to  show  that  each  call  to  Tree-Augment"*"  (and,  analogously,  to  Tree- 
Augment"),  results  in  a  small  number  of  pivots.  More  precisely,  we  have  to  prove  two  claims. 
First,  we  have  to  prove  that  the  intersection  R{T)r\T  is  empty,  then  there  exists  an  arc  {v,w) 
entering  R{T)  such  that  /-^(t;,  w)  <  0.  Then  we  have  to  show  that  there  will  be  only  a  small 
number  of  pivots  (actually,  at  most  m)  needed  to  produce  a  tree  such  that  R{T)  fl  T  is  not 
empty. 

Note,  that  replacing  the  calls  to  Augment  in  the  description  of  the  Fast-Excess-Scaling 
algorithm  (see  Figure  6)  by  calls  to  Tree- Augment"*"  and  Tree- Augment",  causes  Tree- 
Augment"*"  to  be  called  with  parameters  5  =  5/(^A),T  -  Tf(^A). 
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Lemma  4.5  Let  T  be  a  dual  feasible  tree,  f  be  a  pseudoflow  that  is  zero  on  all  arcs  not  in 
the  tree,  s  be  a  node  in  5/(^^A)  and  let  R{T)  denote  the  set  of  nodes  reachable  from  s  over 
residual  arcs  in  the  tree  T.  If  R{T)  Pi  T/(^A)  =  0  then  at  least  one  of  the  tree  arcs  {v,w) 
entering  R{T)  has  f^{v,w)  <  0.  An  analogous  statement  holds  fort  G  Tf{^^A)  instead  of  s. 

Proof:  We  shall  only  prove  the  first  statement.  By  the  definition  of  i2(T)  the  tree  arcs  leaving 
R{T)  are  not  in  the  residual  graph.  This  implies  that  no  arc  in  T  leaves  R{T)  and  also  no  flow 
leaves  or  enters  fl(T).  Consequently, 

veR(T)  v€R(T) 

R{T)  n  T/(iA)  -  0  and  e/(s)  >  ^A  imply  that  Et,6fl(r)  ^/(w)  >  0.  In  particular,  it  follows 
that  R{T)  ^  V.  Now  consider  the  arcs  {vi,wi), .  . .  ,{vk,Wk)  of  the  tree  T  entering  R{T).  The 
sets  V  \  Hr,{v,,w,)  for  f  =  1, .  .  . ,  it  together  with  R{T)  partition  V.  Therefore 

i=l,...,fct,«//T.(„,,u.,)  veR(T) 

The  above  two  equations  and  inequality  X!)t;eR(r)  ^ji^)  >  0  imply  that  J2veV\Hr  (v  w  )  ^^^^  ^  ^ 
for  at  least  one  index  i.  Hence  f    {vi,  Wi)  <  0.       | 


Lemma  4.6  During  the  execution  of  the  procedure  Tree- Augment"*",  no  arc  deleted  from  T 
will  reenter  the  tree.  An  analogous  statement  holds  for  procedure  Tree- Augment". 

Proof:  We  shall  only  prove  the  first  statement.  The  procedure  Tree- Augment"  can  be  treated 
similarly.  Consider  an  execution  of  the  procedure  Tree- Augment"*".  Let  R{T)  denote  the  set 
of  nodes  reachable  from  s  over  residual  arcs  in  the  tree  T.  The  pivot  step  deletes  an  arc  (v,  w) 
for  some  w  £  R{T)  and  adds  to  the  tree  an  arc  {v',w')  for  some  w'  ^  R{T).  This  implies  that 
the  set  R{T)  is  non-decreasing  during  the  execution  of  the  procedure,  and  also  implies  that  an 
arc  which  was  deleted  from  the  tree  can  not  be  added  back.       I 

Lemmas  4.5  and  4.6,  and  Theorem  4.4  imply  the  following  Theorem. 

Theorem  4.7  The  uncapacitated  transshipment  problem  can  be  solved  by  the  dual  network 
simplex  algorithm  in  a  sequence  o/0(nm  log(min{n,  5}))  simplex  pivot  steps. 

Proof:  Lemma  4.5  implies  that  the  simplex  algorithm  presented  above  will  continue  until  it 
finds  an  optimal  transshipment.  Lemma  4.6  implies  that  an  augmentation  is  done  (in  the 
accompanying  Fast-Excess-Scaling  algorithm)  after  a  sequence  of  at  most  m  pivot  steps. 
This  and  Theorem  4.4  imply  that  the  number  of  pivot  steps  throughout  the  algorithm  is 
C)(mnlog(mm{fl,n})).       | 
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A  straightforward  implementation  of  the  algorithm  takes  0{Tn)  time  per  pivot  step  (this  is 
how  long  it  takes  to  choose  the  minimum  reduced  cost  arc  leaving  the  set  H-j  iy^\).  We  get  an 
0(nm^  log(min{fl,  n}))  overall  running  time.  Notice,  that  this  running  time  is  worse  then  the 
one  proved  in  Theorem  3. 16  unless  m  is  close  to  n. 

Using  the  equivalence  of  the  minimum-cost  flow  and  the  uncapacitated  transshipment  prob- 
lems mentioned  in  Section  2  we  get  the  following  corollary. 

Corollary  4.8    The   (capacitated)  transshipment  problem  can   be  solved  by  the  dual  network 
simplex  algorithm  in  (^(m^  log  min{S,  n})  pivot  steps. 
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